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EXPERIENCES WITH TWO-EQUATION TURBULENCE MODELS 


Ashok K. Singhal, Yong G. Lai, and Ram K. Awa N95- 27894 

CFD Research Corporation 
Huntsville, Alabama 


OUTLINE 


• Introduction to CFDRC 

• Experiences with 2-Equation Models 

Models Used 
Numerical Difficulties 
Validation and Applications 
- Strengths & Weaknesses 

• Answers to Three Questions (Posed by Workshop 
Organizing Committee) 

1 . What Are Your Customers Telling You? 

2. What Are You Doing In-House? 

3. How Can NASA-CMOTT Help? 


INTRODUCTION TO CFDRC 


• Young and Energetic (Turbulent) Organization, Dedicated to the 
Continuous Process of Advancement and Effective Transfer of CFD 
Technology 

- PRO IFCT^ 

• TWO TYPES OF ACTIVITIES: 
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INTRODUCTION TO CFDRC (Continued) 


• Objective User of Turbulence Models 

(0, 1 , and 2 Equation Models, RSM and LES) 

• Humble Developer, e.g. Monte Carlo Joint Scalar PDF 

• Active Participant in Recent Small Eddies of Turbulence, e.g. 

- Stanford Endeavor: "Collaborative Testing of 
Turbulence Models" 1989-1993 

National Workshops at: NASA MSFC, LeRC/CMOTT, 
etc. 1987-1994 

- ASME/Fluids Engineering Division, Biathlon, Lake 
Tahoe, June 1994 


TWO-EQUATION MODELS USED 

• Standard k-e Model (Launder & Spalding, 1974) 

• Low-Re k-e Model (Chien, 1982) 

• Extended k-e Model (Chen & Kim, 1987) 

• Multiscale k-e Model (Kim & Chen, 1988) 

• RNG-Based k-e Model (Yakhot et. al. 1993) 

• 2-Layer k-e Model (Rodi, 1991) 

• k-e 4 " 4 " Models 

• k-co Model (Wilcox, 1991) 

++ Models with Corrections for: Curvature, Rotation, Buoyancy, 
Compressibility, etc. 
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NUMERICAL DIFFICULTIES 


• Positivity of k & e (or co) Is Not Guaranteed in Iterative 
Algorithms 

• Strong Nonlinearity of Source Terms and Coupling Causes 
Numerical Difficulties 

• Inappropriate Specifications of e (or co) at Boundaries or 
in Initial Conditions May Also Cause Divergence 

• Non-orthogonaitiy of Grids Adds to Difficulties 

• Non-smooth Change Over for Two-Layer Model Hinders 
Convergence 


VALIDATIONS PERFORMED 

• Channel and Pipe Flows 

• Backward-Facing Step 

• Turnaround Duct 

• Swirl-Flow Combustor 

• Rotating Disk Cavities 

• Boundary Layers 

• Jets, Wakes, and Mixing Layers 

• Periodic Wakes Behind Bluff Bodies 

Examples of Successes and Failures 

1) Flow Around a Square Cylinder; 2) 180° Square Duct; 3) S-Shaped 
Annular Diffuser; 4) Dump Combustor; 5) Backward Facing Step 
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FLOW AROUND A SQUARE CYLINDER 
Strouhal Number 


Strouhal Number = fH 

u o 

f = Frequency of Vortex Shedding 
H = Obstacle Height 
Uq = Freestream Velocity 


Notes: 


1. Experiments By Durao, Heitor, and Pereira (1988) 

2. Computations with CFD-ACE 

Inlet: 78H Upstream; Outlet: 22H Downstream 

Grid: 120 x 80 

Time Steps: Over 70 Per Time Period 

Ref.: Awa, R.K., Singhal, A.K., Lai, Y.G., “Numerical Simulation Of Periodic 
and 3-Dimensional, Turbulent Rows With CFD-ACE," ASME Ruid Dynamics 
Conference, Lake Tahoe, NV, June 19-23, 1994. 


Model/Expt. 

Time 

Period 

Strouhal 

Number 

Expt. 

7.25 

0.138 

Standard k-e 

7.1 

0.141 

2-Layer k-e 

7.1 

0.141 

RNG k-e 

7.6 

0.132 


FLOW AROUND A SQUARE CYLINDER 
instantaneous Streamlines 


Mid-Cycle 



End of Cycle 
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FLOW IN A 180° SQUARE DUCT 


Computational Domain Static Pressure Along Duct Walls 



Ref.: Awa, R.K., Singhal, A.K., Lai, Y.G., "Numerical Simulation Of Periodic 
and 3-Dimensional, Turbulent Flows With CFD-ACE," ASME Fluid Dynamics 
Conference, Lake Tahoe, NV, June 19-23, 1994. 


FLOW IN A 18tf SQUARE DUCT 
Mean Axial Velocity at 0 = 3 
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S-SHAPED ANNULAR DIFFUSER 


O :D*U of SU«n & Fry(1973) 

: Reynolds Stress Model 

:K-£ Model 




* k-e Model and RNG Model Failed to Predict the Correct Location 
of the Maximum Velocity Downstream 

* Computations with CFD-ACE; Publication Under Preparation 


Confined Swirling Flow 
for a Dump Combustor 


O a : Data of Nejad etal.fi 989) 
: Reynolds Stress Model 



* K-e model failed to preserve the vortex core strength 
near center (see x/h=10 & 18) 


* Computational results to be presented at 
1994 ASME Winter Annual Meeting (Chicago) 
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BACKWARD-FACING STEP 



















EXAMPLE APPLICATIONS 


• Gas Turbine Combustors 

• Liquid Rocket Engines 

• Seals and Bearing Cavities 

• Impellers, Inducers, and Fans 

• 1C Engines 

• CFD Reactors 

• External Aerodynamic Flows 

• Plus Many More 


STRENGTHS & WEAKNESSES 

Strengths of 2-Eauation Models 

• Numerically Economical 

• Easy to Modify 

• Reasonable Applicability Within Engineering Accuracy 
Weaknesses 

• Use of Wall Functions Requires First Grid Outside the 
Viscous Sublayer. This is Difficult to achieve, a Priori 

• Low-Re Approach Does Not Offer Overall Advantage. 

• Two-Layer Approach Needs More Work (e.g. Smoothing) 

• Reynolds Analogy Inadequate for Heat-Transfer 
Applications. 

• Effect of Surface Roughness on Turbulence. 
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CMOTT/CP QUESTIONS 


1 . What Are Your Customers Telling You? 


2. What Are You Doing In-House? 


3. How Can NASA-CMOTT Help? 


WHAT ARE CUSTOMERS TELLING? 


• PLEASE Don't Confuse Us, 

with Additional Models and False Hopes 

• Conclusions (Confusion) Over Last 15- Years 

Use k~e Model, with Wall Functions 

Wall Functions, Oh No!, Never!! 

Use Low-Re k~e,.: Which One?, How?? (Good Questions) 

k-e Is No Gooc(,‘ Neglects Non-lsotropicity, etc., etc. 

Jump on RSM Wagon, Now! 

It Can Take You Anywhere, Eventually!! 

Look How Great is this k~e ++ 

When and How to Use it? (Good Questions) 

Look How Accurate is this Scheme, No Numerical Diffusion. 
Don't Contaminate the Solutions with Turbulence 
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WHAT IS CFDRC DOING? 


• Using What is Available, in Best Possible Ways 


• Listening to Both Sides (Model Developers and Users) 


• Trying to Resist Peer Pressures 


• Struggling to Find Resources for Mundane Goals Such as 
Developing Guidelines for Correct Use of Turbulence 
Models 


HOW CAN CMOTT HELP? 


• CMOTT Has Been Providing Commendable Service in the 
Very Difficult Subject: Turbulence 

• "Turbulence Subprogram" Should Help Further 

• Additional Effort is Needed in Many Areas, Such As: 

Near Wall Treatment 

Effect of Surface Roughness 
Economical Heat Transfer Model 

Documentation of Experiences in: 

a) Model Robustness(ln Addition to Accuracy) 

b) Model Sensitivity to Grid Distribution and Boundary 
Conditions 

- Transition Model (if Possible Suitable for k~e Framework) 
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HOW CAN CMOTT HELP? (Continued 


• NASA-CMOTT Is One of the Few Groups Sustaining 
Momentum for Turbulence Modeling. 

• It Is In Unique (Privileged) Position for Embracing the 
Challenge of Developing Specific Recommendations 
(Guidelines) For: 

a) Selection of Adequate Models for Different Class 
of Problems 

b) Correct Use of Each Model 

• The Task Is Difficult But Practical 

• Select Fewer Roads, Post Milestones, and Go Further 

• Move An Inch Closer to Users 


CMOTT 
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